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The structural, vibrational, energetic and electronic properties of hydrogen at the stoichiometric 
RuO2(110) termination are studied using density functional theory. The oxide surface is found to 
stabilize both molecular and dissociated H2. The most stable configuration in form of hydroxyl 
groups (monohydrides) at the undercoordinated O br surface anions is at low temperatures accessed 
via a molecular state at the undercoordinated Ru cus atoms (dihydrogen) and a second precursor in 
form of a water-like species (dihydride) at the O br sites. This complex picture of the low-temperature 
dissociation kinetics of H2 at RuO2(110) is in agreement with existing data from high-resolution 
energy-loss spectroscopy and temperature programmed desorption. Hydrogen adsorption at O br 
sites increases the reactivity of the neighboring Ru cus sites, which are believed to be the active sites 
in catalytic oxidation reactions. 
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I. INTRODUCTION 

The widespread use of transition metal (TM) oxides, 
for example in applications for catalysis, electrochem- 
istry, gas sensors, and corrosion/ wear protection is an 
increasing source of motivation for fundamental research 
on this material class. An important goal in such stud- 
ies, focusing on the surface functionality of oxides, is to 
establish atomic-scale insight into their surface structure 
and composition, as well as their interaction with gas 
phase speciesi*^. Despite notable efforts secure knowl- 
edge is still rather scarce. This holds even for well- 
defined single-crystal surfaces under the controlled con- 
ditions of ultra-high vacuum (UHV), mostly due to the 
structural complexity of oxides and to their often insu- 
lating nature which hampers the use of electron spec- 
troscopy techniques 3 . 

With respect to these issues crystalline RuC>2 repre- 
sents a rather nice choice for a suitable benchmark model 
system. Not only is it one of the few metallic TM oxides 
and its rutile bulk structure of modest complexity*, but 
with e.g. a reported high catalytic activity in oxidation 
reactions^Sii and being discussed as playing a sensitive 
role in Pt-Ru based direct methanol fuel cells^ it is also 
sufficiently interesting from an applied perspective. Es- 
pecially the low energy RuC>2(110) surface has recently 
received particular attention 6 : 7 : 9 : 10 ! 11 : 12 ! 13 : 14 ! 15 ! 16 . Focus- 
ing mostly on the CO oxidation reaction quite some de- 
tailed understanding (experimental as well as theoreti- 
cal) on the fundamental interaction of O and CO with 
this model surface has emerged from these studies. From 
the potential interest for both catalytic and fuel cell ap- 
plications it now appears natural to extend this knowl- 
edge also to the interaction of hydrogen with RuO2(110). 
From a chemical point of view this ubiquitous gas phase 
species is expected to form strong bonds particularly with 
the oxygen anions at the oxide surface^. Intentionally or 
unnoticed hydrogen could therefore be present as sur- 
face species in a wide range of conditions, significantly 
influencing the functionality in the targeted application. 



In fact a noticeable effect of hydrogen contamination on 
the CO turnover numbers has been discussed in a recent 
experimental studyi 7 -. 

On a microscopic level Wang et al. have pro- 
vided detailed kinetic and vibrational data on the 
low-temperature hydrogen adsorption at RuO2(110) in 
UHV 18 . Using high-resolution electron energy- loss spec- 
troscopy (HREELS) and temperature programmed des- 
orption (TPD) they identified both a molecular and a dis- 
sociated hydrogen state, the latter exhibiting vibrational 
properties of a water-like species with, however, a pecu- 
liarly blue-shifted scissor mode. Motivated by these spe- 
cific results we set out to systematically investigate the 
structural, vibrational, energetic and electronic proper- 
ties of hydrogen at the stoichiometric RuO2(110) surface 
using density-functional theory (DFT) 19 . The stability 
of molecular and dissociated hydrogen is first discussed 
at either of the two prominent adsorption sites exhibited 
by this oxide surface (Section IVA and IVB). Then higher 
coverages involving simultaneous occupation of both sites 
are addressed (Section IVC). The detailed picture on the 
low-temperature dissociation kinetics of H2 that emerges 
from the synergetic interplay between computations and 
experiments is intimately connected with the interesting 
property of RuO2(H0) to simultaneously sustain both 
molecular and dissociated hydrogen states, as has already 
been described briefly in a preceding communication 18 . 



II. COMPUTATIONAL DETAILS 

The DFT calculations were performed within the Full- 
Potential Linear Augmented Plane Wave (FP-LAPW) 
me ^-j 10( j2rj 4 2i i 22 usm g ^he generalized gradient approx- 
imation (GGA) 23 for the exchange-correlation func- 
tional. The FP-LAPW basis-set parameters are as fol- 
lows: R^ U T =1.8 bohr, R^ /T =1.0 bohr, R^ /T =0.6 bohr, 
E™{ x =20.25 Ry, E^* x =400 Ry, wave function expansion 
inside the muffin tins up to l^/ ax =12, and potential ex- 
pansion up to lmax = &- The Brillouin zone integration 
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employed a (4 x 9 x 1) Monkhorst-Pack grid with 18 k- 
points in the irreducible Brillouin-zone (IBZ) for (1 x 1) 
surface unit-cells, and a (4 x 4 x 1) grid with 8 k-points 
in the IBZ for (1x2) cells. 

We stress that the very short O-H bonds represent a 
formidable challenge to electronic structure theory cal- 
culations. With respect to the FP-LAPW method em- 
ployed this translates into the necessity to use rather 
small non-touching muffin-tin spheres. Then, the con- 
vergence behavior with respect to the interstitial plane 
wave cutoff is rather slow, requiring to use much higher 
E™4e than the - 17 - 19 Ry typical for late TM oxide 
studies involving muffin-tin sphere sizes only dictated by 
the O-metal bonds. Correspondingly, we tested this con- 
vergence by increasing the plane wave cutoff up to 36 
Ry (!), as well as by employing denser k-meshes up to 
36 k-points in the (1 x 1) IBZ. From these detailed tests 
(cf. the appendix), we conclude that the structural and 
vibrational properties of the systems addressed in the 
present study are well converged within 0.02 A and 8% 
respectively at the above stated E™£ x =20.25 Ry, which 
correspondingly was chosen as our standard cutoff. Rel- 
ative binding energy differences between configurations 
involving the same number of O and H atoms (e.g. when 
comparing tilted and untilted geometries) are similarly 
well converged to within 0.05 eV at this cutoff. For ab- 
solute binding energies on the other hand, in particu- 
lar when they involve strong O-H bonds like in hydroxyl 
groups, the chosen, already rather high cutoff is how- 
ever still not enough. Where such numbers were required 
within a 0.1 eV/H atom numerical accuracy for our phys- 
ical argument and are then listed in this manuscript, we 
correspondingly ran subsequent calculations at a very 
(say for routine calculations still prohibitively) high cut- 
off of E™{,j,=30 Ry employing the structure relaxed be- 
fore at E^-^20.25 Ry. Further, we also checked on the 
uncertainty introduced by the use of the approximate 
exchange-correlation functional by performing a number 
of calculations employing also the local-density approx- 
imation (LDA)24 and will comment on the differences 
between results obtained within LDA and GGA below. 

The RuO2(110) surface was modelled by a three tri- 
layer 0-(Ru2C>2)-0 periodic slab as detailed befor o 12 ! 13 , 
using a vacuum region of about 13 A to decouple the 
interactions between neighboring slabs in the supercell 
geometry. All structures were fully relaxed by a damped 
Newton scheme until the residual forces on the atoms 
were less than 50meV/A, keeping only the atomic posi- 
tions in the central trilayer at their fixed bulk positions. 
Test calculations employing 5 and 7 trilayer slabs and 
relaxing also deeper lying layers showed no significant 
structural changes beyond the topmost trilayer, neither 
was there an influence on the atomic surface geometries 
as obtained with the three trilayer calculations. At this 
point we would further like to emphasize that the struc- 
tural relaxation allowed for any symmetry breaking at 
the surface. This was found to be crucial to obtain the 
correct energetics and structures, which often involve sig- 
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FIG. 1: Computed GGA energy contours as a function of 
Ru02 bulk lattice parameters a and c (with optimized internal 
parameter u). The energy zero corresponds to the optimized 
values of a = 4.52 A and c = 3.13 A marked with '+' and 
higher contour levels are labelled in the figure in meV per 
Ru02 formula unit. Additionally shown are the bulk lattice 
parameters as determined by experiments ( '*' V |27|. 

03 )■ The bulk lattice parameters obtained by an earlier 
DFT-GGA pseudopotential study are a — 4.65 A and c = 
3.23 A 9], outside of the range shown in the figure. 

nificant tilting of the surface groups. 

For the calculations of the vibrational modes of the 
various surface species, the dynamical matrix was set up 
by displacing each of the involved surface atoms from 
their equilibrium positions in 0.04 A steps. Anticipating 
a good decoupling of the vibrational modes due to the 
large mass difference between Ru and O/H, the positions 
of all atoms in the substrate below the adsorption site 
were kept fixed in these calculations. The normal modes 
were then obtained by subsequent diagonalization of the 
dynamic matrix. 



III. BULK RU0 2 , CLEAN (110) SURFACE, 
FREE H 2 AND H a O 

Summarizing and extending the results of our earlier 
publications on RuOzISiiiiiiiSiili we first briefly describe 
the obtained properties of the bulk and the clean (110) 
surface, as well as free H 2 and H 2 molecules, as far 
as they are relevant for the understanding of the hydro- 
genated surface discussed below. 

Ru0 2 crystallizes in the rutile structure, described by 
lattice parameters o and c, as well as one internal de- 
gree of freedom u specifying the positions of the O an- 
ions within the bulk unit-cell^. In order to determine 
the bulk equilibrium lattice constants we scanned a grid 
in a and c near the experimental values with 0.02 A and 
0.015 A steps respectively at a high interstitial plane wave 
cutoff of 30 Ry. At each grid point the internal pa- 
rameter u was further optimized minimizing the com- 
puted forces on the O atoms in the unit-cell. The result- 
ing energy landscape is shown in Fig. ^ while Table ^ 
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TABLE I: R11O2 lattice constants (a, c and u) and bulk mod- 
ulus (-Bo) as determined within LDA and GGA, full-potential 
(FP) and pseudopotential (PP) calculations. Experimental 
values are from X-ray diffraction (XRD) and low- energy elec- 
tron diffraction (LEED) experiments. 





a (A) 


C (A) 


u 


Bo (GPa) 


This work (FP-LDA) 


4.42 


3.05 


0.306 


352 


This work (FP-GGA) 


4.52 


3.13 


0.306 


294 


[25] (PP-LDA) 


4.56 


3.16 


0.307 


283 


[9] (PP-GGA) 


4.65 


3.23 


0.305 




[26] (XRD) 


4.492 


3.106 


0.306 


270 


[27] (XRD) 


4.49f 


3.106 






[28] (LEED) 


4.5f 


3.18 






[9] (LEED) 


4.51 


3.23 








FIG. 2: Side view of the stoichiometric RuO2(110) surface ter- 
mination explaining the location of the two prominent adsorp- 
tion sites corresponding to under-coordinated surface atoms: 
bridging oxygen O br and coordinatively unsaturated (cus) 
ruthenium Ru cus (Ru = light, large spheres; O = dark, small 
spheres) . 



lists the determined bulk lattice parameters within the 
LDA and GGA. Comparing with available bulk diffrac- 
tion studies2SiS£ we obtain very good agreement within 
the GGA. Interestingly, two earlier DFT pesudopotential 
calculations (LDA 25 and GGA 9 .) yield larger lattice pa- 
rameters that deviate by more than 0.1 A from our LDA 
and GGA values, respectively. We particularly checked 
on the GGA value by Kim et al, but find their geom- 
etry to be about 0.4 eV per Ru02 formula unit higher 
in energy than our ground state geometry. Tentatively, 
we take this 0.1 A and 0.4 eV as a rough estimate of the 
structural and energetic error introduced by the pseu- 
dopotential approximation when describing late transi- 
tion metal (TM) oxide compounds. Finally, we notice 
that the c lattice parameter measured in two independent 
low-energy electron diffraction (LEED) experiments 9 ^ 
deviates from the bulk diffraction results by about 0.1 A. 
In both LEED experiments thin RuO2(110) films grown 
on Ru(0001) were investigated, where the c lattice pa- 
rameter corresponds to the size of the surface-unit cell in 
the [001] direction. Apparently, a thin film (with an in- 
commensurable structure to that of the metal substrate) 
prefers a slightly different geometry than a bulk crystal. 



Depending on the oxygen-content in the surrounding 
gas phase, either a stoichiometric or an oxygen-rich ter- 
mination is stabilized at the RuO2(110) surfac o 12 i 16 . Fig- 
ure [21 shows the surface geometry of the stoichiometric 
termination exhibiting two kinds of surface species, the 
nearest neighbor shell of which has been reduced by the 
creation of the surface: A twofold coordinated bridging 
oxygen O br (threefold coordinated in the bulk) and a 
fivefold coordinated ruthenium atom Ru cus (sixfold co- 
ordinated in the bulk). The only other surface species 
present, 3f , still maintains its bulk- like threefold coor- 
dination to in-plane Ru atoms. The oxygen-rich termi- 
nation differs from this geometry only by extra oxygen 
atoms adsorbed on-top of the Ru cus atoms. In the follow- 
ing we will focus exclusively on the interaction of hydro- 
gen with the stoichiometric termination, attempting to 
make contact with the existing data from ultra-high vac- 
uum (UHV) experiments 18 , where the stoichiometric ter- 
mination is the standard surface produced after a high- 
temperature anneal to 600 KSuLiSiii. The discussion of 
hydrogen interaction with the oxygen-rich termination is 
deferred to a consecutive publication. 

Due to the cutting of bonds the surface atoms re- 
lax. Most prominently we find that the O br atoms move 
inwards, reducing their bond length to the underlying 
Ru atoms (henceforth denoted Ru br ) to 1.91 A (bulk: 
1.99 A), and, correspondingly, the inter-layer distance is 
reduced to 1.09A (12% smaller than bulk). These find- 
ings are in very good agreement with a recent LEED 
study which determined a bond length of 1.94 A and 
a first layer contraction of -13% at RuO^llO) 3 -, and 
they are similar to the geometric changes observed at 
the isostructural TiO2(110) surface^ 9 .. Although not as 
pronounced also the under-coordinated Ru cus atoms re- 
lax a bit inwards, thereby inducing a buckling within 
the first trilayer plane of 0.16 A (LEED: 0.18 A). Con- 
trary to other transition metal oxide surfaces like e.g. 
AI2O3 (OOOl)^ 9 . the structural relaxations are therefore 
rather small, and in particular damp away rapidly: No 
significant deviations from the bulk-like positions are 
found for atoms below the first trilayer, neither in LEED, 
nor in our thicker slab test calculations. 

Finally, we summarize in Table ILTI the computed bind- 
ing energies, bond lengths and frequencies for gas phase 
H2 and H2O. For both molecular and atomic calcula- 
tions we employed the same muffin-tin spheres as de- 
tailed before for the slab calculations, and allowed for 
non-spherical densities by reducing the symmetry. The 
total energies for the isolated, spin-polarized atoms are 
obtained by adding to the total energy value from a non- 
spin-polarized FP-LAPW calculation a constant spin- 
polarization energy of 1.52 eV (O) and 1.10 eV (H) taken 
from a relativistic atomic DFT calculation 20 . With re- 
spect to structural and vibrational properties we obtain 
very good agreement with the experimental data, as well 
as with previous DFT studies^S*^. The slow conver- 
gence of the O-H bond energetics already described in the 
preceding section exists similarly for the free molecules, 
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TABLE II: Structure parameters (bond length d and angle), 
binding energy Eb, and vibrational frequencies v for gas phase 
H2 and H2O. Compared are the computed values for LDA and 
GGA (at two different LAPW cutoffs) with the corresponding 
experimental data (with zero-point energy removed). 
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189 
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which is why we list the computed values both at the rou- 
tine 20.25 Ry and at the very high LAPW cutoff of 30 Ry 
(used lateron to obtain quantitative binding energies as 
discussed above) . With respect to the latter cutoff we as- 
sess the numerical convergence of the gas phase binding 
energies to be within 0.2 eV. This points then at a siz- 
able error (i.e. overbinding) compared to the experimen- 
tal values when oxygen is involved even within the GGA, 
which is a well known result^. Although some error can- 
cellation occurs in the computation of binding energies 
at surfaces (as illustrated in the appendix), we conclude 
that a cautious reasoning is necessary when judging on 
the endo- or exothermicity of adsorption with respect to 
the gas phase molecules. To this extent most of our phys- 
ical arguments will instead rather be based on binding 
energy differences, which are more accurately described 
(±0.05 eV, cf. the appendix). 

IV. THE HYDROGEN ATED RUO 2 (110) 
SURFACE 

In the following sections we systematically discuss the 
hydrogen interaction with the various surface species 
present at the stoichiometric RuO2(H0) termination. 
The energetics will be described with respect to the afore- 
described free H2 molecule, where a negative binding 
energy denotes endothcrmicity with respect to the gas 
phase species and the bare surface, i.e. a metastable situ- 
ation. As we find the interaction with the threefold coor- 
dinated in-plane 3f to be energetically very unfavorable 
(even the formation of a monohydride is endothermal by 
« — 0.3eV/H atom), the discussion will concentrate on 
the two under-coordinated surface species, i.e., Ru cus and 



O br . At first we discuss hydrogen first to be present at 
one of these sites only (Section IVA and IVB). And then 
we consider simultaneous adsorption at both sites (Sec- 
tion IVC). Consistent with our previous publications we 
will employ a short hand notation to characterize the 
manyfold of studied geometries, indicating first the oc- 
cupancy of the bridge site and then of the cus site, e.g. 
(OH) br /H2 US for a configuration with an OH-group at the 
bridge site and a H2 molecule at the cus site. 

Lateral interactions between functional groups at di- 
rectly neighboring cus and bridge sites are implicitly con- 
tained within our calculations employing (lxl) surface 
unit-cells. With these cells only integer multiples of one 
monolayer (ML) hydrogen coverage can be studied (1 ML 
defined as monoatomic occupation of all sites of one type 
(br or cus)). Further reaching lateral interactions, e.g. 
towards a moiety at the same site type either along [110] 
(at a distance of 6.4 A) or along [001] (at a distance of 
3.1 A, see Fig. |3J) could lead to the formation of more 
dilute superstructures with fractional MLs hydrogen cov- 
erage. From a systematic study of oxygen adsorption at 
RuO2(110) we find such lateral interactions to be rather 
small at this rather open oxide surface^ and don't expect 
this to be significantly different for hydrogen, in partic- 
ular along the longer [110] direction of the surface unit- 
cell, where the next-nearest site would be 6.4 A away. 
Correspondingly, we only test for lateral interactions in 
selected configurations employing larger (1 x 2) cells, al- 
lowing us to model structures where then only every sec- 
ond site is occupied along the [001] direction (see e.g. 
Fig. □ below). 



A. Hydrogen at Ru cus 

We start with the molecular adsorption of H2 at the 
under-coordinated Ru cus site (O br /H™ s ). This analy- 
sis was motivated by the aforementioned recent UHV 
HREELS experiments that attributed a weak peak at 
367 meV to the stretch mode of a molecular hydrogen 
species at the surface with a TPD activation energy of 
about 0.3 eV— . Relaxing a H2 molecule in (1 x 1) cells 
from a height at about 2 A atop the cus sites our cal- 
culations indeed find such a (meta)stable species with a 
computed binding energy of +0.32 eV/H 2 with respect 
to molecular H2. The resulting adsorption geometry is 
shown in Fig. |3| together with the calculated vibrational 
modes. At a height of 1.85 A the H2 molecule lies parallel 
to the surface above the cus sites (side-on configuration) . 
Interestingly, we find almost no corrugation of the poten- 
tial energy surface with respect to an azimuthal rotation 
of the flat-lying H2 molecule: The optimal bond orien- 
tation about 30° degree from the [110] direction is only 
by insignificant 3 meV more stable than any other orien- 
tation, i.e. the H2 behaves essentially like a free-rotator 
(helicopter mode), as also reflected by the very low ro- 
tational vibration frequency of 12meV, cf. Fig. [3J The 
similarly low in-plane translational modes further point 
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stretch 4 1 3 in eV rot ati en 1 68 m eV rotation 1 2 si eV 




trans. 93 raeV trans, 5 6 meV trans 38 raeV 

FIG. 3: a) Top and side view of the adsorption geometry of 
1 ML H2 at Ru cus (Ru = light, large spheres; O = dark, small 
spheres; H = white, small spheres) . Additionally shown is the 
size of the (1 x 1) surface unit-cell. Note that the azimuthal 
direction of the plotted H2 is arbitrary, as it is freely rotating, 
b) Six vibrational modes of the adsorbed H2 molecule. The 
directions and lengths of the arrows indicate approximately 
the directions and relative magnitudes of the displacements 
of the atoms. 



at a rather expressed derealization of the H2 molecule 
parallel to the surface, particularly in the [110] direc- 
tion, i.e. approaching the neighboring bridging oxygens. 
Despite this the H2 bond length is with 0.81 A notice- 
ably stretched compared to the free gas phase molecule 
(0.75 A). This is consistent with the rather high binding 
energy and the significantly decreased stretch vibration 
(413 meV compared to 538 meV in the gas phase). 

These properties reflecting a moderately strong inter- 
action are very similar to molecular H2 at late transition 
metal surfaces (like e.g. Pd(100)2£) if adsorption is re- 
stricted to take place at the on-top site. Of course, at 
the latter surfaces the on-top site is only a local "con- 
strained" minimum and H2 would dissociate towards 
higher-coordinated hollow sites. The different geometry 
of RuO2(110) doesn't offer such sites, thus stabilizing 
the molecular adsorption at Ru cus . In this respect we at- 
tribute the non-dissociative interaction of H2 with these 
sites more to a geometry effect, rather than to an elec- 
tronic effect, i.e. compared with the weak physisorption 
of H2 at noble metal surfaces like Ag(lll)^. 

The relatively strong interaction (electron polarization 
and bond formation) is also obvious in the "difference 
density" plot shown in Fig. The "difference density" is 
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0.05 




FIG. 4: Difference density plot for H2 adsorbed at Ru cus . 
The contour plot depicts the plane parallel to the H2 molec- 
ular axis and normal to the surface. Areas of electron ac- 
cumulation and depletion have positive and negative signs 
respectively, contour lines are drawn at 0.01 e/A 3 intervals. 



obtained by subtracting from the electron density of the 
adsorbate system that of the clean surface and that of the 
free molecule (the latter two with the same interatomic 
distances as found for the adsorpbate system)^ . A rather 
strong polarization of the flat-lying H2 is apparent, with 
electron density accumulation on the substrate side of the 
molecule and depletion on the vacuum side. Due to the 
interaction with the H2 the back-bond of the Ru cus atom 
to the underlying O substrate atom is slightly weakened, 
the Ru atom moves 0.03 A upwards and thereby reduces 
the clean-surface-buckling in the topmost trilayer. 

The structure of the difference density plot suggests 
the major interaction in the occupied states to be due 
to a hybridization of Ru-d z 2 and liVcr orbitals. Analyz- 
ing the computed local density of states we indeed find 
the bonding governed by this and a smaller hybridization 
of the Hu-d xz /d yz with the H2-CT* orbitals. The bond- 
ing of the H2 molecule to the cus site can therefore be 
understood within the familiar donation/back-donation 
picture^, where the hybridization with the H2-CT orbital 
causes a small H2— >Ru charge transfer that is counter- 
acted by some back-donation of electronic charge from 
the metal to the H2-CT* level, strengthening the coupling 
to the substrate while weakening and elongating the H-H 
bond. 

Contrary to the situation at most more reactive TM 
surfaces this back-donation is however not strong enough 
to completely dissociate the H2 molecule. This is supple- 
mented by the surprising result that atomic hydrogen 
is not stable at the Ru cus site: We compute only an 
endothermal binding energy of — 0.33eV/H atom with 
respect to 1/2 H2 at 1 ML H-coverage, i.e. in (1 x 1) 
unit-cells. Checking whether H bonding might become 
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more favorable in more dilute superstructures we also em- 
ployed (1 x 2) cells to model a 0.5 ML H-coverage with H 
only at every second site along the [001] direction. With 
-0.20 eV/H atom with respect to 1/2 H2 the binding en- 
ergy is still endothermal and only slightly changed, re- 
flecting the small lateral interactions at the RuO2(110) 
surface. Expecting no further changes in binding energy 
for even more dilute H-phases we therefore conclude that 
only molecular hydrogen may be stabilized at the cus 
sites. 

Forming the basis for e.g. relations between hetero- 
geneous and homogeneous catalysis it is interesting to 
compare these findings for the Ru atom at the surface 
of an oxide with the hydrogen bonding to TM atoms 
in other frameworks like e.g. at the surface of met- 
als or in a TM complex. As already mentioned, cou- 
pling of hydrogen to TM surfaces is generally associated 
with the dissociation of the ligancUi. The observation 
of non-dissociative chemisorption of H2 has so far been 
restricted to a few exceptional case a 42 i 43 , mostly con- 
nected with a prior saturation of the most reactive sites 
at the surface with atomic hydrogen. Concomitantly, we 
compute e.g. the bonding of atomic H at Ru(0001) to 
be exothermic at least up to 1 ML coverage. Molecular 
precursors have more been identified at noble surfaces 
like Ag(lll)iS, yet then physisorbed and certainly not 
exhibiting such a strong activation of the H2 bond as 
expressed at RuO2(110) with the significant bond elon- 
gation and downshift of the stretch frequency. These 
findings resemble much more the data from organometal- 
lic complexes: For so-called r) 2 -H2 (dihydrogen) sin- 
gle metal atom complexes, in which the H-H bond re- 
mains intact 4 *, TM-H2 bond energies in the range of 0.1- 
0.3eV/H2 are estimated 4 -, significant red-shifts of the 
stretch frequency 4 ^* 4 ^ and bond elongation up to 0.9 A 
are reported^ 4 ^. In fact, neutron scattering experiments 
furthermore indicate rapid rotation of f? 2 -H2 ligands with 
an activation energy of less than about 10meV— , just as 
we find for the molecular hydrogen at Ru cus . Even the 
donation/back-donation bonding model is analogously 
discussed for the TM complexes 4 ^. Yet, there (just like 
at TM surfaces) the back-donation may also break the 
H-H bond, and often H2- as well as H-ligands are at- 
tached to the same metal center and can even exhibit 
continuous changes between both configurations 4 ^* 4 ^. At 
RuO2(110), in contrast, the bonding to the surround- 
ing oxide apparently depletes the electron density at the 
under-coordinated Ru cus atom already in such a way, 
that the back-bonding only activates the H2 bond, but 
may no longer break it. 
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FIG. 5: a) Top and side view of the adsorption geometry 
of 1 ML H at O br , computed in (1 x 1) unit-cells, b) Six 
vibrational modes of the formed hydroxyl group. 
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B. Hydrogen at O br 

For H2 molecules approaching the RuC>2 surface at the 
under-coordinated O br atoms we find that they will ei- 
ther slide towards the Ru cus sites or that they are repelled 
into the gas phase. In no case (i.e., testing many H2 ori- 



FIG. 6: Difference density plot of a hydroxyl group at the 
bridge sites (OH) br /— . The contour plot depicts the plane 
parallel to [110] and normal to the surface. Areas of electron 
accumulation and depletion have positive and negative signs 
respectively, contour lines are drawn at 0.01 e/A 3 intervals. 
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FIG. 7: a) Top and side view of the adsorption geometry of 
1 ML H at O br , computed in (1x2) unit-cells with alternating 
tilting, b) Top and side view of the adsorption geometry of 
0.5 ML H at O br . 



entations, as well as lower coverages in a (1 x 2) cell) 
did. we observe spontaneous dissociation above the O br 
site. This leads us to conclude that molecular adsorption 
primarily takes place over the aforediscussed cus sites. 
Atomic hydrogen, on the other hand, binds at the O br 
sites, i.e., forming a surface hydroxyl group. The opti- 
mized geometry of the 1 ML (OH) br /— phase computed 
in (1 x 1) cells is shown in Fig. |SJ The strong binding 
within the OH-group weakens the bond to the underly- 
ing Ru br substrate atoms and elongates it significantly 
from 1.91 A to 2.06 A. The surface OH- group itself has 
a bond length of 0.97 A, nearly identical to the O-H dis- 
tance in H2O, and is inclined towards the [110] direction 
with a 40° angle with respect to the surface normal. The 
computed binding energy is +0.89 eV/H atom stronger 
than in molecular H2, and the energy gain by the tilting 
is 0.1 eV/H atom compared to the higher-symmetry, up- 
right configuration. Again, both the computed stretch 
frequency, as well as the binding energy are found to be in 
good agreement with the recent HREELS and TPD data 
from the hydroxylated surface^, and the strong binding 
is also nicely visible in the difference density plot shown 
in Fig. El 

With only a 0.1 eV energy difference between tilted and 
upright position, the hydroxyl groups will at finite tem- 
peratures frequently swing from one orientation to the 
other. Checking whether this may occur in concerted 
wave-like motions along a chain of bridge sites, we re- 
peated the calculation with 1 ML coverage, but now in 
a (1 x 2) cell with each OH-group alternatingly tilted 
in one or the other direction as shown in Fig. \7^,. The 
only structural difference obtained is a somewhat larger 
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FIG. 8: a) Top and side view of the adsorption geometry of 
0.5 ML water-like species at the bridge sites, oriented along 
the [110] direction, b) Corresponding nine vibrational modes. 
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FIG. 9: Top and side view of the adsorption geometry of 
0.5 ML water-like species at the bridge sites, oriented along 
the [001] direction. 



tilt angle of about 60° degree, while bond length and 
binding energy remain to within 0.01 eV/H atom vir- 
tually unchanged. Testing for further reaching lateral 
interactions by decreasing the total coverage to 0.5 ML 
as shown in Fig. \7]p, we again find the binding energy 
within 0.02 eV/H atom degenerate to the previous two 
cases, leading us to conclude that each OH-group tilts 
and swings essentially independent of the others, and of 
the total coverage which may thus easily reach the full 
1 ML with each O br atom hydroxylated. 
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The low temperature HREELS and TPD experiments 
by Wang et ali& show that the coverage of dissociated 
hydrogen can be increased above 1ML. There is evidence 
for a hydrogen state at O br , that exhibits frequencies 
resembling a water-like H-O-H configuration, i.e. with 
two H atoms per bridge site (dihydride). Surprisingly, 
the measured scissor mode for this functional group is 
with 231 meV higher than the one of gas phase water 
(exp: 198 me V, calc: 189 meV, see Table ITT)l . whereas 
intuitively a red-shift would be expected: Considering 
that the molecule-surface interaction weakens the intra- 
molecular bonds and normally widens the H-O-H angle, 
both these effects would rather tend to make the bend- 
ing mode softest. Trying to address this puzzling find- 
ing we computed such water-like configurations at the 
bridge sites first at 1ML H20-coverage in (1 x 1) cells 
and with the water-axis oriented once along [001] and 
once along [110]. Both configurations exhibit only a very 
low stability and can thus not account for the experi- 
mental water-like species with a desorption temperature 
around 300 K 1 -. We therefore proceeded to relax the 
same two configurations at lower coverages in (1 x 2) cells 
with water-like species now only at every second bridge 
site as shown in Figs. |S] and While the model with 
orientation along the [001] axis is still slightly endother- 
mal, the binding energy of the other orientation shown 
in Fig. [S] finally turns out already at least exothermal 
by +0.48eV/H 2 . Still, as almost expected the calculated 
scissor mode for this configuration (just as much as the 
one of the other three models) is with 168 meV signifi- 
cantly lower than the one of free water, and thus in strong 
disagreement with the experimental data. 

Recalling the strong tilting of the mono-hydride group 
at O br shown in Fig. [5]we then tested to similarly tilt the 
whole water-like species of the last most favorable model 
and ended up with the geometry shown in Fig. 1101 that 
is apparently separated from the up-right configuration 
by a sizable energy barrier (i.e. neither configuration re- 
laxes automatically into the other one). The tilt not only 
significantly increases the binding energy of this configu- 
ration by 0.08 eV to +0.56eV/H 2 , but also the computed 
scissor mode is with 211 meV 12% higher than the one 
of a free water molecule (189 meV, see Table lH|) . Also 
the other computed vibrational modes listed in Fig. 110b 
are now in reassuring agreement with the experimental 
HREELS data (exp: stretch 436 meV, scissor 231 meV, li- 
bration 110 meV 76meV, translation 59meV 28meV) 18 . 

The question remains why this tilted geometry gives 
rise to the counterintuitively blue-shifted scissor mode. 
Inspecting the relaxed geometry displayed in Fig. ITUI 
in more detail we find the water-like H-O-H bond an- 
gle with 108° slightly increased compared to free water 
(calculation: 103°, see Table lll| . lending more towards 
the argument favoring a weakening of the bending mode. 
However, this angle together with the overall tilting of 
the whole functional group causes the lower OH-bond 
of the water-like species to end up almost parallel to 
the surface (only 8° to the surface plane), bringing the 
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FIG. 10: a) Top and side view of the adsorption geometry of 
0.5 ML water-like species at the bridge sites, oriented asym- 
metrically along the [110] direction, b) Corresponding nine 
vibrational modes. 




FIG. 11: Difference density plot of the asymmetric water- 
like species of Fig. 1101 at the bridge sites. The contour plot 
depicts the plane parallel to [110] and normal to the surface. 
Areas of electron accumulation and depletion have positive 
and negative signs respectively, contour lines are drawn at 
0.01 e/ A 3 intervals. 
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terminal H atom closer to the neighboring Ru cus atom 
(2.72 A compared to 3.19 A in the upright model shown 
in Fig. |HJ. This suggests a possible interaction between 
these two species that is indeed verified by inspecting 
the computed difference density plot shown in Fig. 1111 
Although small, the apparent hydrogen bridge bonding 
possible in this tilted configuration hinders any move- 
ment away from Ru cus parallel to the [110] direction and 
thus naturally stiffens not only the scissor mode, but also 
the corresponding librational mode (calc: 119 meV, exp: 
110 meV, compared to "typical" values around 80meV, 
cf. Figs. EUlandlHJ). 

This intricate coupling between neighboring bridge and 
cus sites fits nicely into the picture of the low temperature 
dissociation kinetics of H2 at RuC>2(110) that emerges 
from the data presented in Sections IVA and IVB, as 
well as from the detailed UHV-HREELS experiments by 
Wang et al. 18 . Although the thermodynamic ground 
state for hydrogen at this stoichiometric surface is given 
by the strongly bound hydroxyl-groups at the bridge sites 
(binding energy +0.89eV/Hatom)A&, direct dissociation 
over these sites is apparently inhibited at least at low 
temperatures by an energy barrier. Thus, we tend to 
conclude on the following scenario: When H2 interacts 
with RuC>2(110), it is at first bound in molecular form 
at the cus sites (binding energy +0.16 eV/H atom). Al- 
though the H2 bond is weakened in this process, the cus 
sites can not induce its complete cleavage. This is instead 
achieved via the water-like metastable configuration at 
the bridge sites (binding energy +0.28 eV/H atom) that 
may be accessed from the molecular H 2 state at the cus 
sites (as indicated by the low translational mode in [110] 
direction, cf . Fig. [3J . From there the hydroxyl groups are 
finally formed, presumably by activated hydrogen diffu- 
sion. This interpretation of H™ 55 as a necessary precursor 
state to dissociation is also supported by new low tem- 
perature HREELS and TDS experiments which report a 
suppressed population of the water-like species at bridge 
if the cus sites are first blocked by CO molecules^. Only 
if the water- like species are allowed to form, the hydroxyl 
groups result from moderate heating to 350 K 18 . 



C. Hydrogen at both sites: Higher coverages 

Having discussed the lower coverage adsorption up to 

1 H2-ML for both sites separately, we now proceed to 
higher total coverages involving hydrogen at both sites. 
Also with various hydrogen functional groups present at 
the bridge sites we still find atomic hydrogen at cus to 
be always unstable. Correspondingly we restrict our de- 
tailed discussion to molecular H2 at cus and different 
hydrogen populations at the bridge sites. Starting with 
the water-like species at bridge, the (1 x 1) configura- 
tion (H 2 0) br /ri2 US corresponding to a total coverage of 

2 H2-ML is computed to be endothermal, reflecting pre- 
sumably already an oversaturation of the surface with 
hydrogen. This improves in more dilute (1 x 2) con- 




FIG. 12: Top and side view of the adsorption geometry of the 
3H-ML (OH) br /Hl us phase. 



figurations with (H20) br and occupying only ev- 

ery other site. The average binding energy for both a 
checkerboard arrangement and with both species occu- 
pying directly neighboring sites turns then out exother- 
mic by about +0.2eV/H 2 , with a slight preference for 
the checkerboard arrangement in which the two hydro- 
gen species maximize their mutual distance. This average 
binding energy is, however, still lower than the average 
binding energy we obtained for the (lxl) O bl /H™ 8 phase 
(+0.32eV/H 2 , cf. Section IVa) which also corresponds to 
a total coverage of 1 H2-ML. 

This leaves as interesting higher coverage phase be- 
yond 1 H2-ML only the possibility to combine H™ 13 with 
hydroxyl groups at the bridge sites. For this remain- 
ing combination even a dense (lxl) arrangement of 
(OH) br /H| us at 1.5H 2 -ML total coverage turns out very 
stable with an average binding energy of +0.6eV/H2. 
This points at the possibility that after the aforedis- 
cussed formation of hydroxyl groups at bridge also the 
molecular state at the cus sites could simultaneously be 
populated upon continued hydrogen uptake. In other 
words that the RuO2(110) surface offers the fascinating 
property that hydrogen may coexist both in the dissoci- 
ated monohydride and in the non-dissociated dihydrogen 
state. Inspecting the corresponding geometry shown in 
Fig. El the first striking effect of the simultaneous oc- 
cupation of bridge and cus states is that the pronounced 
tilting of the hydroxyl group, cf. Fig. [SJ has disappeared. 
Next, the H2 at the cus sites is apparently more activated 
compared to the situation discussed in Section IVA when 
the bridging oxygen atoms were bare: The bond length is 
increased from 0.81 A to 0.83 A, and the molecule resides 
at 1.80 A height, i.e. 0.05 A closer to the surface. 

Interestingly, the formation of the hydroxyl group at 
the bridge sites seems to influence the bonding properties 
at the neighboring cus sites, rendering the latter some- 
what more reactive. The mechanism with which this hap- 
pens is nicely identified on the basis of the difference den- 
sity plot shown in Fig. Plotted are the induced elec- 
tron density variations arising in the (OH) br /H2 US phase 
when the hydroxyl group is formed. Apart from the obvi- 
ous significant density rearrangement at the bridge sites 
themselves (which is very similar to the one shown in 
Fig. El for the lower coverage (OH) br /— phase), also 
some variations can be observed at the distant cus sites: 
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FIG. 13: Difference density plot of the high coverage (1 x 
1) (OH) b 7H£ us phase shown in Fig. [H Displayed are the 
induced electron density changes in the system when hydrogen 
forms a hydroxyl group at the bridge sites. The contour plot 
depicts the plane normal to the surface that cuts through a 
bridge site and its next-nearest neighbor cus site (this way 
passing very close a 3f atom and the positions which are 
nearest to the 3f are indicated by arrows). Areas of electron 
accumulation and depletion have positive and negative signs, 
respectively, contour lines are drawn at 0.01 e/A 3 intervals. 



electron density is depleted around the Ru cus atoms and 
clearly populates the H2 cr*-orbital, thereby weakening 
the molecular bond. At the same time the H2-metal 
bond is strengthened, so that the absolute binding en- 
ergy is only slightly altered (0.23eV/H2 compared to the 
0.32eV/H2 without the presence of the hydroxyl group, 
cf. Section IVA). Still, the molecule is more activated 
and adsorbs closer to the surface. 

We interpret the surface's capability to give a higher 
back-donation to the cus bond upon population of hy- 
droxyl groups at bridge as arising from an interesting 
bond order propagation effect. Due to the newly formed 
hydroxyl bond the O br atoms become less bound to 
the underlying Ru br atoms, as reflected by the consid- 
erable increase in bond length from 1.91 A to 2.01 A. 
Having thus lost a bit of their optimum bonding envi- 
ronment, the Ru br seek to fortify their remaining back- 
bonds, among others also to the directly coordinated in- 
plane 3f atoms, cf. Fig. [3 The resulting electron 
density rearrangement at the latter atoms can nicely be 
discerned in Fig. 1131 which displays a plane connecting 
next-nearest neighboring bridge and cus sites and thus 
passes also closely by the 3f sites. With this slightly 
strengthened bond to the Ru br atoms the 3f atoms in 
turn adapt by weakening their bonds to the Ru cus atoms. 
Due to this propagation of bond order the latter atoms 
find themselves in a neighborhood that is less electron 
demanding, allowing for an increased back-donation into 
the H 2 -Ru cus bond. 

While we find further reaching lateral interactions at 
this surface (e.g. the aforediscussed low coupling between 
hydroxyl groups at neighboring cus sites) to be rather 
small, the bonding properties at nearest-neighboring 



bridge and cus sites are thus obviously under a non- 
negligible mutual influence. This suggests that one could 
attempt to tune the reactivity of either site by a con- 
trolled population of the respective other site. Particu- 
larly for oxidation reactions at RuC>2(110) several stud- 
ies have already emphasized the key role played by the 
cus s itesfii2iiLiiiiiii. Correspondingly, modifications at 
the bridge sites like the present decoration with hydroxyl 
groups could have a noticeable impact on the overall cat- 
alytic activity of the surface - not only because of a 
possible site-blocking of bridge sites, but also because 
of an intricate tuning of the electronic structure at the 
cus sites. So far, all reaction mechanisms discussed at 
this RuO2(H0) surface have been found to be initial- 
state dominated, i.e. the reaction barriers scaled with the 
binding energies of the adsorbed reactantsii 3 -. The above 
described increased electron-donation ability of the cus 
sites upon hydroxylation of the bridge sites might there- 
fore likely lead to a different reactivity of this oxide sur- 
face. Concomitantly we notice that a promoting effect of 
small amounts of hydrogen on the CO oxidation reaction 
over polycrystalline Ru0 2 has already been reported in 
a recent experimental study, proposing this material as a 
suitable candidate for the technologically important low- 
temperature CO oxidation in humid air—. The present 
work identifies an increased reactivity of the cus sites af- 
ter hydroxylation of the bridge sites at RuO 2 (110). Fur- 
ther work involving atom-specific Surface Science exper- 
iments on the reactivity of hydrogenated Ru0 2 surfaces 
is now required to check if this can be exploited not only 
for total oxidation reactions, but possibly even more im- 
portant for partial oxidation, where a controlled tuning 
of the reactive sites might be crucial to obtain a high 
selectivity. 



V. SUMMARY 

In conclusion we presented a detailed ab initio study of 
the energetic, electronic, structural and vibrational prop- 
erties of hydrogen at the stoichiometric RuO2(110) ter- 
mination. A different interaction with the two underco- 
ordinated, prominent adsorption sites at this surface is 
found: At Ru cus only a molecular state can be stabilized, 
while the thermodynamic ground state is represented by 
hydroxyl groups involving the O br surface atoms. The 
results strongly suggest that the low-temperature H2 dis- 
sociation takes place via a precursor where molecular 
H2 is bound at the Ru cus site. From there a water- 
like dihydride state at O br can be accessed that leads to 
the final hydroxyl groups. This emerging picture of the 
low-temperature dissociation kinetics has been developed 
in a synergetic interplay with the experimental group 
of Jacobi and Ertl, and is in full agreement with their 
HREELS and TDS datai&, in particular with the pecu- 
liar blue-shifted scissor mode measured for the water-like 
species at O br . 

Upon further hydrogen adsorption both the molecular 
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H2 state at the Ru cus site and the dissociated mono- 
hydride state at bridge-site surface oxygen atoms be- 
come populated. The formation of the hydroxyl groups 
is hereby found to intricately influence the reactivity at 
the neighboring cus sites, allowing for an increased back- 
donation further activating the H™ 3 bond. This modifi- 
cation of the bonding properties at cus by hydrogen deco- 
ration at bridge is attributed to a bond order propagation 
mechanism, possibly special to this metallic oxide. It is 
argued that the resulting possibility of fine tuning the cus 
site reactivity by controlled modification of the bridge 
site population could be of relevance for catalytic appli- 
cations, in particular partial oxidation reactions where a 
precise tuning of the bond strengths could be crucial to 
obtain high selectivities. 
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Appendix 

In order to ensure the needed numerical accuracy we 
particularly performed calculations with larger intersti- 
tital plane wave cutoffs and denser k-point meshes - the 
latter two being the most influential parameters in the 
FP-LAPW basis set. Some of these tests are sketched in 
this appendix. Specifically, we increased the cutoffs up 
to 36 Ry in several steps and doubled the k-point mesh 
to 36 k-points in the IBZ for (1 x 1) surface unit-cells. 
Table Hill lists the resulting bond lengths, bond angles, 
and binding energies for three characteristic phases in- 
volved in the present study: the bare RuO2(110) sur- 
face (O br /— ), a hydroxyl-group at bridge ((OH) br /-) 
and molecular hydrogen at cus (O hl /H^)- It can be 
seen that the structural parameters of all three phases 
are well converged already at the lowest cutoff listed 
(20.25 Ry), let alone that the employed k-mesh has any 
influence. However, the absolute binding energies are not 
yet fully converged at this low cutoff, primarily due to the 
need to use very small muffin-tin spheres as discussed in 
the text. From the convergence behavior along the se- 
quence 20.25 Ry-^25Ry^30.25Ry^36Ry we conclude 



that only at a high cutoff of 30.25 Ry the latter quanti- 
ties are converged to within 0.1 eV/H atom. Fortunately 
enough, relative energetic differences between similar ge- 
ometries involving an equal number of O and H atoms 
converge much more rapidly. This is illustrated by a 
comparison of the binding energy of the hydroxyl group 
at bridge either in an upright or in the tilted configura- 
tion (listed in Table IIII|) : at all tested cutoffs between 
20.25 Ry and 36 Ry this energetic difference is O.leV/H 
atom, constant to within 0.01 eV/H atom. Correspond- 

TABLE III: Computed bond lengths (d in A), tilt angle (in 
°) and absolute binding energies (E b in eV) as a function of 
interstitial plane wave cutoff and number of k-points for three 
phases characteristic for the present study. 
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E b (0) 


20.25 Ry, 18 k 1.91 








2.54 


25 Ry, 18 k 1.90 








2.46 


25 Ry, 36 k 1.90 








2.43 


30.25 Ry, 18 k 1.91 








2.43 


36 Ry, 18 k 1.91 








2.37 
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-Ru br 




tilt 


E b {U) 


20.25 Ry, 18 k 2.06 




0.97 


40 


0.71 


25 Ry, 18 k 2.08 




0.98 


45 


0.84 


25 Ry, 36 k 2.08 




0.98 


46 


0.82 


30.25 Ry, 18 k 2.07 




0.98 


46 


0.89 


36 Ry, 18 k 2.07 




0.97 


46 


0.91 


br /H cu S dQbr 


-R.u br 


dH-Ru cus 




E b (U 2 ) 


20.25 Ry, 18 k 1.92 




1.89 


0.81 


0.38 


25 Ry, 18 k 1.92 




1.89 


0.81 


0.35 


25 Ry, 36 k 1.92 




1.89 


0.81 


0.34 


30.25 Ry, 18 k 1.92 




1.89 


0.81 


0.32 


36 Ry, 18 k 1.92 




1.89 


0.81 


0.30 



ingly, we employed the manageable cutoff of 20.25 Ry 
for structural relaxations, vibrational calculations and 
when judging on the relative energetic sequence of simi- 
lar structures. Only when absolute binding energies con- 
verged to within 0.1 eV/H atom are required, did we run 
additional calculations at 30 Ry, but with fixed geometry. 
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